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Abstract 

Using a minimal model based on the continuum theory of a 2D hard-core/square-shoulder 
ensemble, we reinterpret the main features of cluster mesophases formed by colloids with soft 
shoulder-like repulsive interactions. We rederive the lattice spacing, the binding energy and 
the phase diagram. We also extend the clustering criterion [Likos, C. N., et al. Phys. Rev. E, 
2001, 63, 031206; Glaser, M. A., et al. EPL 2007, 78, 46004] to include the effect of the hard 
cores, which precludes the formation of clusters at small densities. 

Introduction 

Genetically, neutral colloidal particles attract each other via van der Waals and Casimir forces, 
and measures must be taken to prevent aggregation of clusters in both experimental systems and 
technological applications. One might be led to conclude that purely attractive hard-core particles 
could only form close-packed lattices or glassy messes. But the phenomenon of aggregation in 
colloids is not restricted to particles that attract each other. Over the past two decades clustering 
in purely repulsive pair potentials has been explored in some detail to find that it is distinguished 
by emerging order not seen in attractive particles. In particular, it has been established that at large 
enough density, clumping repulsive colloids can form superstructures with large voids on the or- 
der of many particle diameters. 1 Indeed, Malescio and Pellicane 2 demonstrated that even simple 
hard-core, square-shoulder potentials led to clustering purely on energetic grounds. The cluster 
morphologies and the clustering criterion itself have been studied theoretically using a range of 
approaches including liquid-state theory, 3 lattice theory, 4,5 density functional theory, 6 and contin- 
uum models, 5,7 and the predictions of the different approaches are remarkably consistent. Equally 
unequivocal are the results of numerical studies, mostly using Monte Carlo methods 5,8 and direct 
search of minimal-energy configurations using genetic algorithms.- - — 

It is worthwhile at this juncture to step back and develop a minimal set of simple rules to ex- 
pose the mechanism of lattice formation in these systems. It is our hope that these rules will make 
clear the essential ingredients needed for pattern formation in this large class of purely repellent 
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systems. Moreover, we develop these ideas in the context of real-space potentials and interactions 
yet recapitulate the clustering criterion of the more technical (though precise) Fourier-based anal- 
yses.-^' 5 Designing potentials in real space offers a more intuitive route to rational self-assembly 
and connects directly with, for instance, laser-trap and depletion based potentials. 

In this note, we do this by using the continuum T = model of the stripe phase formed by 
particles interacting via the hard-core/square-shoulder pair potential 



U(r) 



°°, r < O 

£, O <r <A (1) 
0, r> A 



where o and A > o are the diameters of the core and the shoulder, respectively, and £ is the 
shoulder height. 12,13 Characterized by a one-dimensional density modulation, the stripe phase is 
mathematically the most transparent of all cluster morphologies and our real-space analysis may 
be easier to visualize than reciprocal- space arguments. 3 " 5 Moreover, by focusing on energy rather 
than on free energy we emphasize that the entropy does not promote clustering. These results 
are then extended to explore the generic phase diagram of the cluster-forming system at finite 
temperatures. 



Clusters at T = 

In the following, we first evaluate the energy of a two-dimensional hard-core/square- shoulder sys- 
tem with given hard-core diameter o and shoulder diameter A at a fixed average number density, 
and we minimize it with respect to intra-cluster density, cluster size, and lattice spacing. We con- 
sider the simplest cluster morphology, the one-dimensional square-wave of uniformly populated 
parallel stripes of width d and particle-free gaps of width £ — d; the lattice spacing is t. In the 
continuum model suitable in the large shoulder-to-core limit where A is sufficiently larger than a, 
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the energy per particle can be related to the average overlap area per particle defined as 



1 

co = — 



2d Jo 




(2) 



where 0(r) is the Heaviside step function describing the shape of the shoulder potential and ri 
and T2 are the locations of particles 1 and 2, respectively, measured from the origin at an edge of 
the stripe in question. The first integration over y\, the distance of particle 1 from a stripe edge, is 
over the stripe containing particle 1 and the integration over Y2 goes over all stripes, co captures the 
interaction of the shoulders of particles, whereas the hard-core repulsion is treated in a mean-field 
approximation by demanding that the number density within stripes p st ripes be no larger than 



corresponding to the close-packed hexagonal arrangement of the particles' hard-disk cores. This 
approximation is applicable in the broad shoulder regime A>(J and at densities large enough that 
the stripes are sufficiently wider than the core diameter so that speaking of "intra-cluster packing" 
of particles has meaning. 

In terms of CO, the average energy per particle reads 



But the number density of particles within stripes depends on their width relative to lattice spacing 
d/£, which represents the fraction of the total area that is occupied by the stripes. In terms of the 
average density p, 
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(3) 



Pcp = 



V3o 2 



E — £Pstripes&'- 



(4) 




(5) 



Since the energy is to be minimized at fixed p rather than at fixed pstripes it is convenient to intro- 
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duce the scaled average overlap area 



CO 



which includes all dependence of 

E = epQ. (7) 

on d and £. For the hard-core/square- shoulder pair interaction, £1 can be computed analytically but 
the result is too cumbersome to be of interest here. 
We will now establish some rules of thumb. 



I. Clusters are close-packed In [figure] [1][]1 we plot the reduced scaled average overlap area 
Cl = Q./X 2 as a function of reduced lattice spacing £ = £/X for several stripe widths dji. The 
cluster-free, unmodulated phase corresponds to £ = d = and its nature depends on density: As 
argued below, its phase sequence includes the expanded fluid, the expanded hexagonal crystal, the 
condensed fluid, and the condensed hexagonal crystal. But since we treat the hard-core part of 
the pair interaction in a mean-field fashion, worrying only about the average number of neighbors 
within the reach of a particle's shoulder, the exact nature of the unmodulated phase is not crucial 
to understand why clustering takes place. 

The uniform, unmodulated phase can be interpreted as a stripe morphology with a very fine 
density modulation such that the lattice spacing and the stripe width are much smaller than the two 
characteristic length scales of the pair potential, X and o. In this limit, Cl = n/2 irrespective ofd/£ 
which tells us that for d,£ <^ X the average energy per particle is E = nepX 2 /2. As £ is increased, 
Cl oscillates around it/ 2, reaches a global minimum at £ « 1.2, and then grows monotonically to 
saturate at a value of n£/2d. The large /behavior is a signature of macroscopic phase separation: 
At fixed d/£, states with reduced lattice spacing £ beyond ~ 1 correspond to thick stripe widths d 
with an ever smaller number of particles residing at the boundary of the stripes. The average energy 
per particle is gradually dominated by that of the particles well within the bulk of the stripes, which 
reads E = 7T£p str i pes /l 2 /2 = K£pX 2 £/2d because the density within the stripes is larger than the 
average density by a factor of £/d. This result gives Cl(£ — > °°) = n£/2d. 



Q 




Figure 1 : Reduced scaled average overlap area as a function of the reduced lattice constant i = i/X 
for d/i = 0.1(0.1)0.9; for the sake of clarity, only curves corresponding to d/i = 0.1,0.5, and 0.9 
are labeled. The global minimum of the reduced scaled overlap area Q. at £~ 1.2 is deepest for 
vanishingly small d/£, which suggests that the equilibrium stripes are as compact as possible, i.e., 
close-packed. The square density waves representing stripes with dji = 0.1,0.5, and 0.9 (top) 
schematically depict the stripe morphologies at £ = 0.5, 1, 1.5, and 2; the second column where 
£ = 1 corresponds to a lattice spacing exactly equal to the shoulder width. 



The most important feature of this diagram is that the depth of the minimum decreases with 
increasing relative stripe width d/i. This means that at any given average density p = p s tripes^/^> 
the system will select the state with the smallest possible d/i at the expense of the density within 
stripes. For example, the curves corresponding to d/i = 0.1 and 0.2 represent two possible states 



of the system of a fixed average density p. According to | [figure] [!][]![ the absolute minimum of 
the former is lower than that of the latter, which means that of the two stripe phases in question, 
the d/i = 0.1 state minimizes the total energy E = epX 2 Cl. But since the relative stripe width of 
the d/i = 0.1 state is half of that of the d/i = 0.2, the corresponding density within stripes must 
be twice as large as in the latter state. In other words, the ground state geometry of the stripes 
minimizes d/i, restricted only by the close-packing limit forbidding stripes with p str ipes = pi/d 
beyond the close-packed density p cp . We conclude that at T = the optimal stripes are made of 
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close-packed particles so that p s tri P es = Pep and the reduced average number density 



n=f- (8) 

Pep 



coincides with djl. 



II. Lattice spacing weakly depends on density An additional feature of the graphs in [figure] [ 1 ] [] 1 
is that the equilibrium reduced lattice spacing, £ eq , depends only weakly ond/£. It is apparent that 
£ eq is largest at half-filling where it reaches 1.217, which nicely agrees with the value of 1.223 
predicted by the lattice theory. 5 It can be shown that the dependence of £ eq on n is well described 
by a parabola symmetric about n = 0.5. Even in the infinitely-dilute and close-packed limits, l sq 
tends to 1 — but this result is of no physical consequence. As shown below, the binding energy 
vanishes in these two limits so that the clustering mechanism is not acting and the corresponding 
equilibrium reduced lattice spacings are irrelevant. 

This observation is made even more compelling by employing observation I: Since the ground- 
state stripes are close-packed, we may as well switch back to the average overlap area per particle, 
CO and write the energy as E = £p cp co; CO = co/X 2 can thus be regarded as the reduced energy. In 



[figure] [2][]2| we plot cd as a function of reduced lattice spacing I. It is apparent in this presentation 
that the location of the global minimum of CO is little changed over the range 0.1 < n < 0.9 and 
again we conclude that £ eq is roughly independent of n. 



From [figure] [2][]2 we also extract the reduced binding energy Aco defined as the depth of the 



global minimum of co(n) at £ eq relative to the unmodulated phase at £ = 0. The reduced binding 



energy shown in the inset to [figure] [2] [] 2 is a skewed U-shaped function of reduced average 
density which vanishes for n = and n = 1. In an infinitely dilute system, the energy of the 
unmodulated phase itself tends to and no spatial modulation of the density profile can reduce it 
further so that the binding energy is too. On the other hand, a system of average density close to 
p cp cannot undergo but a very restricted spatial modulation (because the density within the stripes 
should not exceed p cp ) and thus the energy gained upon clustering approaches when p — > p cp . 



to 



1.5 - 
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Figure 2: Reduced energy of the close-packed stripe morphology, 00, as a function of reduced 
lattice spacing I for n = 0.1(0.1)0.9. In the unmodulated phase at £ = 0, 00 = nn/2. The inset 
shows the reduced binding energy Aob = G>(leq) — G)(£ = 0) [as illustrated in the G)(n = 0.5) curve] 
which is a skewed U-shaped function of reduced average density and vanishes at n = and n = 1. 



Although 00 and £2 are closely related, they convey a somewhat different message. From 



[figure][l][]l we learned that the particles within the stripes are close-packed, which enabled us to 



directly relate the average density p to relative stripe width d/£. On the other hand, | [figure] [2] [] 2 
exposes the binding energy of the equilibrium stripe morphology more clearly. But as far as the 
magnitude of the equilibrium lattice spacing is concerned, both quantities are equally telling. 



Phase diagram 

Our two observations can be used to qualitatively outline the phase diagram of the cluster phases 
in the temperature-density plane. To this end, we need to study the difference of the free energies 
of the stripe phase and the unmodulated phase, which consists of an energy term AE = ep cp A 2 A(U 
and of an entropic term. The entropy of the two phases depends on their structure elaborated below 
in a semi-quantitative fashion. 

At absolute zero, one of the hallmark features of the stripe morphology is its compact intra- 
stripe structure where the impenetrable hard cores of the particles are packed together as tightly 
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as possible at any reduced average density n. This state is materialized by the hexagonal lattice. 
Because of the robust, density-independent nature of this behavior, we posit that at finite temper- 
atures the intra-stripe density p str ipes should not depend strongly on the reduced average density 
either (though it must be smaller than the close-packing density p cp ) and that the intra-stripe order 
remains hexagonal. 

The structure of the unmodulated phases of hard-core/square- shoulder particles is more com- 
plicated and despite decades of efforts (see, e.g., Refs.— ~— ), their thermodynamics remains only a 
partly solved problem. For the purpose of present discussion, it suffices to note that at low temper- 
atures where they compete with cluster morphologies, the phase sequence of unmodulated phases 
consists of 4 variants, the 2 low-density and high-density phases being characterized with little 
and sizable overlap of particles' shoulders, respectively. At very small densities, the particles form 
an expanded fluid of disks of shoulder diameter A (schematically shown in [figure] [3] []3[ ). As this 



fluid is compressed, it undergoes a transition to the expanded hexagonal crystal of disks of diameter 
A . The location of the transition can be estimated by rescaling the phase diagram of the hard-disk 
system: p e f- e c ~ 0.792(c/A) 2 p C p. 20 Upon further compression, the expanded hexagonal phase 
remelts to avoid close-packing; since the overlap of shoulders is increasingly less unfavorable at 
elevated temperatures, the phase transition density should decrease with temperature. From this 
transition on, the particles behave essentially as hard spheres of diameter o and the transition to 
the condensed hexagonal phase takes place approximately at p c f- C c ~ 0.792p cp .— 

Understanding the main features of the sequence of unmodulated phases helps us to construct 
semi-qualitatively the entropic part of their free energy. We first note that the pressure of the hard- 
disk hexagonal crystal can be roughly regarded as a continuation of the hard-disk fluid branch 20 
and we approximate the excess entropic free energy per particle in both the fluid and the crystal 
unmodulated phase by 

f an 

(9) 



Fg*{n) = k B T 



ln(l — an) 



I - an 

where a = n/2\/3 ~ 0.907. This prediction is based on the 2D Carnahan-Starling type theory for 
the fluid phase.— ^ Admittedly a rough approximation — it does not distinguish between the fluid 



and the crystalline phase, and the hexatic phase is disregarded altogether — this excess entropic 
free energy provides a simple and adequate model across a broad range of reduced densities con- 
sistent with the scope of this analysis. Its main deficiency is the poor description of the crystalline 
phase (whose excess entropic free energy should diverge at reduced average density of n = 1 and 
not at n = l/a = 1.103). Yet we note that in the phase diagram of hard disks, the condensed 
crystalline phase is stable at n > 0.792— so that the discrepancy is limited only to large densities. 

Using the model F ex [Eq. (??)], we can outline the excess entropic free energy for the un- 
modulated phases. In [figure][3][]3[ we plot the excess free energy per particle of the expanded 



and the condensed phases, F^ p (n) = F e c x ond ((A/a) 2 rc) and F£° nd (n). The former diverges at 
n = (o/X) 2 /a which corresponds to close-packed disks of diameter A. This divergence is, of 
course, unphysical because the expanded hexagonal lattice remelts upon compression so that in 
this regime, the true excess entropic free energy interpolates between F^ p and F^° nd . 
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Figure 3: Excess entropic free energies per particle of the expanded and the condensed unmodu- 
lated phases derived using a Carnahan-Starling type theory (solid lines; the expanded phase cor- 
responds to XJo = 2). The divergence of F^ p at n = 0.276 is unphysical; instead of approach- 
ing the shoulder-to- shoulder close-packed structure, the system melts to form the condensed fluid 
phase. The dashed line represents a qualitatively correct interpolation between the expanded and 
the condensed branch. — The schematics illustrate the structure of the 4 unmodulated phases: The 
expanded fluid, the expanded crystal, the condensed fluid, and the condensed crystal. Full circles 
indicate the hard cores of particles whereas the shaded coronas represent the shoulders. 



The exact shape of the interpolating F ex is not known. But our mean-field model is designed to 
work best for broad shoulders and in this case, the existence of the expanded phases is restricted 
to reduced average densities below (a/A) 2 / a, i.e., to very small n. Being interested in the overall 
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behavior of the system, we may approximate the excess entropic free energy by the condensed 
branch alone. Then the total free energy difference can be constructed from overlap energy differ- 
ence AE(n), the excess entropic free energy of the stripe morphology F e c x ond (n str i pes ), and the excess 
entropic free energy of the unmodulated phase F™ nd (n): 

AF(n,n stripes ,T) = AE(n)+F™ nd (n stlipes )-F™ nA (n) 
= ep C pA 2 < n stripes A<a(n) 



""stripes . OCn 



ln( 1 - awstdpes) - i h ln( 1 - an) 



1 - (XHstripes 1 - an 



(10) 



Here n s tripes — Pstripes/ Pep is the reduced density within stripes, which should not depart much from 
1 and must decrease with temperature, and 

- w o» 



epcpA 2 
is the reduced temperature. 

The clustering criterion The most important features of the total free energy difference are the 
negative skewed U-shaped energy term, whose exact dependence on the reduced average density 



shown in the inset to [figure] [2] [] 2 can be well fitted by 

A£^-0.65ep cp A 2 n(l-n) 2 , (12) 

and the positive excess entropic free energy difference proportional to temperature which mono- 
tonically decreases from a finite value at n = to at n = n str i pes . We approximate the temperature 
dependence of the reduced density within stripes by a linearly decreasing function 

«stripes(0 = 1 -CX. (13) 
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In the following, we choose c = 6: At the largest reduced temperature where the stripe morphology 
is stable, this gives n s tripes ~ 0.9 which is plausible. 



In [figure] [4] []4i, we plot AF for several values of reduced temperature T. As t is increased, 
the reduced average density range where the stripe morphology is stable gradually shrinks and 
at a large enough T, the stripe morphology is disfavored at any n. Thus the phase diagram is 



characterized by a dome-like region of stability of the stripe morphology shown in [figure] [4] [] 4 3. 
Except at very large n where the differences between the unmodulated and the stripe phase are 
increasingly smaller and our model is inaccurate, the shape of the phase boundary reproduces well 
the clustering criterion obtained in terms of the lattice theory 5 which states that instability occurs 
when 

n(l — n) 

— - > const. > 0. (14) 

T 

The agreement is indeed remarkable although our phase diagram lacks the symmetry about the 
half-filling point n = l/2 encoded in Eq. (??). Needless to say, the model can be refined by better 
estimating the entropic part of the excess free energy. 

Extending the clustering criterion The stripe morphology is subject to two consistency con- 
straints. Firstly, adjacent stripes separated by more than a shoulder width do not interact with 
each other because there is no overlap of particles residing within them; without a restoring in- 
terstripe repulsion, such configurations would spontaneously disintegrate into thinner stripes with 
narrower gaps between them. This means that in the mechanically stable stripe state, the width 
of the particle-free gaps between stripes £ — d should be smaller than the shoulder diameter A. A 



close inspection of [figure] [1][]1 shows that the global minimum at reduced lattice spacing I « 1 
is not compromised by this condition at any reduced average density n. Secondly, the stripe width 
must be larger than the hard-core diameter of the particles, d > a; if not, speaking of close-packed 
stripes does not make sense. These two conditions bracket the physically relevant range of lattice 
spacing £ from top and from bottom, respectively. To compare the two bounds, we replace d by n£ 
so that i) the upper bound £ — d<X becomes l{\ — n) < X wherefrom £ < A/( 1 — n) and ii) the 
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Figure 4: Free energy difference of the stripe and the unmodulated phase for Stripes described by 
Eq. (??) and reduced temperature x = 0,0.005,0.01,0.015, and 0.02 (panel a). As temperature 
is increased, the range of reduced average densities n where AF < and the stripe morphology 
is stable (thick sections of curves) becomes increasingly more narrow. Each curve terminates at 
n * ^ "stripes ( T ) sucn that AF(n > n*) < 0; states beyond this point correspond to average density 
very similar to the density within stripes where the predictions of our model are meaningless. 
— Panel b) shows the temperature-density phase diagram of the hard-core/soft-shoulder stripe 
morphology and the unmodulated phase computed using the mean-field continuum model with 
the model n st ripes(f ) [Eq. (??)]. The stripe morphology is stable in the shaded region whose shape 
agrees rather well with the clustering criterion [dashed line; const, in Eq. (??) adjusted to reproduce 
the slope of the phase boundary at small reduced average densities] . 
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lower bound d > a becomes n£ > a and thus £ > a/n. We thus find that a stable phase can only 
exist if 

n>y-^— , (15) 

the lower limit of stability of the stripe morphology. In view of the nature of the continuum model 
used here, this treatment of the hard-core part of the pair potential is expected to be valid for 
core-to-shoulder ratios X/o sufficiently smaller than 1 . 

Thus the effect of the hard-core part of the pair potential is to disfavor clustering at small 
densities, thereby restricting the validity of the criterion [Eq. (??)] to densities beyond a threshold 
determined by the core-to- shoulder ratio. Since the hard-core interaction is athermal, this condition 



should apply at all T as depicted in [figure] [5] []5 



Conclusions 

The complete phase diagram will include the fluid and one or more crystal lattices as the low- and 
the high-density variants of the unmodulated, non-cluster phase. On top of the stripe morphology, 
in two dimensions there also exist the disk and the inverted disk cluster phase. 5 Just like stripes 
are most stable at about n = 0.4 ( [figure] [4] []4[ ) which corresponds to £ « 1.2 and d = 0.48, the 



disks are expected to be bound most tightly at a similar lattice spacing and disk diameter. This 
configuration of disks will cover a smaller fraction of the plane and the corresponding average 
density will be smaller than 0.4. If we assume that the disks are stable in a dome-like region of 
the phase diagram qualitatively similar to that describing the stripes, the disk dome must peak at 
a density smaller than that of the stripe dome. Conversely, the inverted disk morphology should 
prevail at densities larger than 0.4. Moreover, the disk and the stripe phases occur in both the 
liquid and the solid intra-cluster order 5 and so the generic full phase diagram of a cluster-forming 
ensemble should have a multiple-dome structure curbed by the fluid phase at small densities and by 



the crystal phases at large densities ([figure] [5] []5 ). At large temperatures, the fluid-crystal phase 



transition line must be vertical because for T — > °°, the system reduces to hard disks of diameter a 
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which freeze and melt at p/p cp = 0.779 and 0.792, respectively. 20 

temperature 




density 

Figure 5: Generic phase diagram of cluster- forming repelling particles. At low temperatures, the 
phase sequence includes the fluid phase (F), several cluster morphologies (fluid disks — Df, solid 
disks — D s , fluid stripes — Sf, solid stripes — S s , and inverted disks — ID), and one or more 
crystal phases (X). The vertical boundary of the disk morphologies at small densities indicates 
the restriction imposed by the hard-core part of the pair potential [Eq. (??)]. Regions of phase 
coexistence are not shown for clarity. 



[figure] [5] [] 5 reproduces many features of the phase diagram obtained using a more complete 
treatment of the thermodynamics of a hard-core/square- shoulder system 5 and it bears some simi- 
larity to the phase diagram of the hard-core/linear-ramp system. 19 Although the width of the ramp 
studied in Ref .— is too narrow for fully developed cluster phases, the non-close-packed lattices oc- 
curring in the phase diagram are very reminiscent of the cluster morphologies discussed here, and 



the phase diagram itself has roughly the same multidome shape as that in [figure] [5] []5 We expect 
that for the hard-core/square- shoulder potential with small core-to-shoulder ratio, the agreement 
of the numerically obtained phase diagram with our prediction should be even better. 

The ideas presented here capture the main mechanisms of cluster formation in systems of 
classical repelling particles in a way marked by the appeal of real-space description and by the 
analysis of the density-modulated morphologies across the whole range of lattice spacing. Given 
the seemingly counterintuitive behavior of particles with shoulder-type pair interaction, we hope 
that our rederivation will clarify the details of the more elaborate studies. 
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